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Abstract 

We investigate the effect of mesonic fluctuations on the pion propagator in the 
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■ Am the pion properties change only quantitatively we encounter strong instabilities 
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1 Introduction 

It is generally believed that chiral symmetry, which is an approximate symmetry of quan- 
tum chromodynamics, is spontaneously broken in vacuum. Together with the linear sigma 
model the Nambu-Jona-Lasinio (NJL) model plays the role of a prototype for this mech- 
anism. The original papers were published long time before QCD was developed Jl[. By 
the observation that massless pions emerge as a consequence of the spontaneous symmetry 
breaking they had direct impact on the formulation of the Goldstone theorem || . Today 
the NJL-model, which was originally a model of interacting nucleons, is reinterpreted as 
a schematic quark model || |4j]. Being much simpler than QCD, it is an attractive tool 
for studying consequences of chiral symmetry breaking as well as its restoration at large 
temperatures or densities, although its applicability to real processes is limited by several 
shortcomings of the model, in particular the lack of confinement and non-renormalizability. 
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So far, most NJL model calculations have been performed in leading order in 1/N C , 
the inverse number of colors. Some authors have also considered corrections in next-to- 
leading order in 1/N C , e.g. to discuss it — tt scattering |J or meson loop corrections to the 
pion electromagnetic form factor §. With the appropriate choice of model parameters 
chiral symmetry is spontaneously broken in leading order, resulting in the appearance of 
massless pions and a massive sigma meson. If care is taken in the choice of correction 
terms 0, g, the pion emerges massless also in next-to-leading order. However, in a recent 
paper Kleinert and Van den Bossche argue that for the physical value, N c = 3, the 
spontaneous breakdown of chiral symmetry does not occur, as a consequence of strong 
fluctuations, which come about as higher-order corrections in the l/A^-expansion ||. At 
present this paper is subject to many controversial discussions. 

In the present article we want to study the effect of mesonic fluctuations on the pion 
propagator by calculating l/iV c -corrections explicitly. However, instead of trying to give a 
definite answer to the question whether or not chiral symmetry is spontaneously broken, 
we take the point of view that there is no unique solution in a non-renormalizable model. 
At each order one encounters new divergences, which are in principle related to additional 
parameters (e.g. counter terms). In the present case it is rather natural to introduce a 
new cutoff parameter Am to regularize the meson loops || . By varying this cutoff from 
zero to larger values we can smoothly turn on the mesonic fluctuations and study their 
effects. 

Although we cannot directly investigate the question of a possible restoration of chiral 
symmetry within our scheme, one might find hints for an instability. This is a similar sit- 
uation as in many-body systems, where e.g. complex energy eigenvalues in the excitation 
spectrum indicate an instability of the ground state [JT0| . In fact, as we will discuss in 



Sec. 3, at large values of A M the pion propagator receives "unphysical" poles, i.e. poles at 
complex or space-like momenta. We also find negative wave function renormalization con- 
stants corresponding to residues with the wrong sign. The main point we want to make in 
this letter is that the existence of such signals depends on the choice of the meson cutoff 
A M . Obviously, for small values of A M the l/-/V c -corrected results are close to the leading 
order ones and no instability shows up. Ultimately, of course, the cutoff should be fixed 
by a fit to observables, for instance the decay width of vector mesons [ |TT|j . 



2 The Model 

In this section we give a brief outline of our scheme for describing mesons within the NJL 
model in next-to- leading order in 1/N C . To a large extent our model is based on Ref. 0. 



More details will be presented elsewhere ]TT 



We consider the standard NJL Lagrangian for two flavors and three colors with scalar- 
isoscalar and pseudoscalar-isovector interaction: 

C = - m )^ + g [(^) 2 + (^75^) 2 ] • (1) 

Here g is a dimensionful coupling constant of order 1/N C . In the limit m = the 
Lagrangian is symmetric under SU(2) l x SU(2)r transformations. 
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Figure 1: Upper part: 1/N C - correction terms to the inverse meson propagators. The double 
lines denote meson propagators in leading order in 1/N C , which are calculated by iteration 
of quark- antiquark loops as shown in the lower part of the figure. 



In leading order in 1/N C mesons are described by iterated quark- antiquark loops, 
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This is illustrated in the lower part of Fig. |T|. For simplicity we will call Dfj a "prop- 
agator", although strictly speaking it should be interpreted as gM qq D%, where is 
the renormalized meson propagator and gMqq is a meson quark coupling constant. The 
quark- antiquark polarization diagrams are given by 
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with Tm = ilh f° r the pion and = 1 for the sigma channel. In Eq. (||) m denotes the 
constituent quark mass which is a solution of the gap equation 



m = m + 2ig (N f N c ) 



tr 
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(4) 



It is now straight forward to show that for mo = but m/fl ("Nambu-Goldstone mode") 
the right hand side of Eq. (0) has a pole at q 2 = in the pion channel and at q 2 = Am 2 
in the sigma channel, corresponding to = and m a = 2m ]!],(§]. If mo/m is small 
but not exactly zero the pion mass becomes non-zero and is approximately given by the 



Gell-Mann Oakes Renner (GOR) relation |TJ: 

m lfl = -mo{#> 



(5) 



The 1/iVc-correction terms to the inverse meson propagators, are shown in the upper 
part of Fig. [I]. These diagrams contain mesonic fluctuations calculated in leading order 
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in 1/N C , as described above. When iterated together with Il^f the new terms lead to an 
improved meson propagator, 

DMiq2) = i-2 9 (n%tf) + snMtf)) ' (6) 

where 5U(q 2 ) denotes the sum of the correction terms. It can be shown analytically that 
the pion constructed in this way is again massless in the chiral limit |7], |TTJ . Note that we 
do not use the improved meson propagators Dm(q 2 ) for evaluating the correction terms 
5U(q 2 ). Such a selfconsistent procedure spoils the 1/N C counting scheme and leads to 
inconsistencies with chiral symmetry. 

A non-renormalizable model is incomplete without defining how to regularize divergent 
loops. To leading order in l/N c one has to deal with quadratically divergent integrals, the 
qq polarization diagrams (Eq. @) and the gap equation (Eq. (f|)). These divergences can 
be regularized by standard methods, such as the Pauli-Villars scheme or a sharp cutoff 
Q. In next-to- leading order the situation is more involved. As an example we consider 
diagram (a) of Fig. [I]. This diagram contains two intermediate gg-mesons, which have 
to be calculated in a first step. We are then left with three loops: two quark triangles 
and an integration over the four-momentum of the intermediate gg-mesons. It is quite 
natural to regularize the quark triangles in the same way as the qq polarization diagrams 
which enter into the intermediate mesons. In fact, this is necessary in order to preserve 
chiral symmetry. However, there is no stringent reason, why the remaining meson loop 
should also be regularized in this way. We therefore follow Refs. |7| and || and introduce 
a meson cutoff as an independent parameter. 

For computational convenience we work in the rest frame of the meson (i.e. the external 
three-momentum is zero) and regularize the internal meson loops by a 3-dimensional sharp 
cutoff Am in momentum space. On the other hand, since the internal mesons can of course 
not be chosen to be at rest, we should use a covariant scheme for their computation. As 
in Ref. J7| we employ the Pauli-Villars scheme to regularize the quark loops. Because of 
the quadratic divergences we need two regulators. The regulator masses are chosen in the 



standard way as Mi = ^Jm? + and M 2 = ^m 2 + 2A^, with a free parameter A q . This 
"mixed" scheme with a covariant regularization of the quark loops and a non-covariant 
regularization of the meson loops was chosen for entirely practical reasons. However, we 
expect qualitatively similar results for other schemes, e.g. if we take a covariant sharp 
cutoff for both, quark and meson loops. 

In this context we should also comment on Ref. [13] where l/A c -corrections to the 
quark condensate are studied within a NJL-type model with a separable non-local inter- 
action. This interaction generates a form factor at the quark- vert ices and can be chosen 
such that the 3-momentum of each quark is limited to absolute values less than a certain 
cutoff parameter. In this way also the 3-momentum of the meson loops is automatically 
restricted without introducing an additional cutoff. Note that this does not contradict 
the statement that additional regulators are necessary if a non-renormalizable model is 
extended to higher loop orders: Because of the non-local interaction the model of Ref. [|T^] 
is not identical to the NJL model but a modification which contains no longer any di- 
vergences. Although this is a quite appealing feature, the model has the disadvantage of 
being manifestly non-covariant (the form factors depend on the 3- momenta of the quarks). 
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Figure 2: Squared pion mass as a function of the current quark mass mo for different 
meson loop cutoffs: A M = MeV (solid), 500 MeV (dashed), 900 MeV (dashed- dotted) 
and 1300 MeV (dotted). The calculated points are explicitly marked. 



There are similar models with 4-momentum dependent form factors (in Euclidean space) 



14] , |15f but these models have wrong analytical properties, which might be very disturb- 



ing in the context of this article. Therefore - and because of the numerical effort - we 
decided to restrict our investigations on the standard NJL model with the regularization 
scheme described above. 



3 Numerical Results 

We are now ready to study the influence of l/iV c -corrections on the NJL pion propagator 
and related quantities. We begin with the leading order, which corresponds to a meson 
cutoff Am = 0. With A q = 800 MeV, gA 2 q = 2.9 and m = 6.1 MeV we obtain a reasonable 
fit for the pion mass, the pion decay constant and the quark condensate: = 140 MeV, 
f( 0) = 93.5 MeV and (^) (0) = -2 (241.1 MeV) 3 . Here and in the following the superscript 
(0) is used to denote quantities which are calculated in leading order in l/N c . The above 
parameters correspond to a relatively small constituent quark mass of 260 MeV. 

Now we turn on the mesonic fluctuations by taking a non-zero meson cutoff Am- Fig. |] 
displays the squared pion mass as a function of the current quark mass mo for different 
values of A^. Obviously all points which correspond to the same meson cutoff lie almost on 
a straight line through the point (mo = 0, m\ = 0). The latter was calculated analytically 
whereas all other points are numerical results. This demonstrates the consistency of our 
scheme with chiral symmetry and the stability of the numerics. 

Fig. |3|(a) displays the behavior of m%, f% and the quark condensate with increasing 



5 




500 1000 1500 4 8 12 



A M [MeV] q 2 /(140MeV) 2 

Figure 3: (a) The ratios m^/m 2 }^ (solid), f%/f% (dashed), (^pip) / (ipip)^ (dashed- 
dotted) and the combination —mo^^/m 2 / 2 (dotted) as a function of the meson loop 
cutoff Km- (b) Inverse pion propagator, D~ x , multiplied by —2g as a function of the 
^-momentum squared for various meson cutoffs: Km = (solid), 900 MeV (dashed), 
1300 MeV (dashed- dotted) and 1500 MeV (dotted). 

Km- All other parameters are kept constant at the values given above. As one can see the 
mesonic fluctuations lead to a reduction of f n while is increased. There are two 1/N C - 
correction terms to the quark condensate, which have been discussed, e.g., in Ref. fill . 
We find that the absolute value of the quark condensate decreases at smaller values of 
K M but goes up again for A A/ > 900 MeV. 

The leading-order quantities , f 2 ^ and are in almost perfect agreement 

with the GOR relation, Eq. (^). From the (almost) linear dependence of m 2 on m , 
as shown in Fig. |2], one might expect that this is also the case for Km > 0. If one 
carefully expands both sides of Eq. (|5|) up to next-to-leading order in 1/N C , one can show 
analytically that our model is consistent with the GOR relation [pjj . However, for the 
quantities m 2 and f 2 as they follow from the l/iV c -corrected inverse pion propagator, the 
l.h.s. of Eq. (|^) also receives higher-order terms (i.e. terms beyond next-to-leading order) 
in 1/N C , which are not present on the r.h.s. and therefore violate the GOR relation. This 
can be seen in Fig. |3](a), where the ratio of the right hand side and the left hand side of 
Eq. @ is displayed by the dotted line. For K M < 900 MeV the relation holds within 30%. 
However, when the meson cutoff is further increased the deviation grows rapidly. This 
indicates that higher-order corrections in 1/N C become important and our perturbative 
scheme should no longer be trusted in this regime. 

The behavior of the pion mass becomes more clear from Fig. |3](b) where the inverse 
pion propagator D~ l (q 2 ) is plotted as a function of q 2 for different values of Km- In the 
chiral limit all lines in this plot would go through —2gD~ x = at q 2 = 0. For mo > 
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the q 2 = 0- value of the leading-order result (solid line) is shifted up to m /m, the ratio 
of current and constituent quark mass, and the curve crosses the zero-line at a positive 
value of q 2 , corresponding to a finite pion mass. For Ajv/ ^ 900 MeV the main effect of the 
l/-/V c -corrections is a further increase of the q 2 =0- value, while the slope does not change 
significantly (dashed line). This leads to a moderate increase of the pion mass. If we 
continue to increase Am, the g 2 =0-value of —2gD~ l moves down again. However, at the 
same time the function becomes more flat and this causes the accelerated increase of the 
pion mass one sees in Fig. |3|(a). For Am > 1250 MeV the curve even turns around at larger 
q 2 and crosses the zero line a second time (dashed-dotted line). At Am ~ 1350 MeV the 
two poles of D n merge and upon further increasing Am disappear from the positive real 
g 2 -axis (dotted line). 

Obviously the pion becomes unstable if the mesonic fluctuations are too strong. In fact, 
already the second pole of the propagator in the regime 1250 MeV < A A/ < 1350 MeV 
is unphysical because the residue has the wrong sign. This would correspond to an 
imaginary pion-quark coupling constant and to a negative value for f 2 . We found that 
this behavior is mainly due to diagram (b) in Fig. [TJ This diagram has a peak in the 
qq continuum which has an imaginary part with the "wrong" sign. When Am exceeds a 
certain value, this contribution becomes larger than the sum of all other contributions, 
such that the imaginary part of the inverse pion propagator, lmD~ l , gets the "wrong" sign. 
Via dispersion relations this peak can be related to the turn-around of Re/)" 1 below the 
qq threshold which is responsible for the second pole. Because of the unphysical features 
of the second pole we used the first one to determine the pion mass for Am = 1300 MeV 
in Fig. 0. 

It is quite reasonable that the instabilities of the pion propagator indicate an instability 
of the underlying ground state with a spontaneously broken chiral symmetry against 
mesonic fluctuations. On the other hand it is clear that our model is not applicable to 
study the process of chiral symmetry restoration itself. Since 1 /Ac-corrections have been 
built in only perturbatively, the range of validity of the model is restricted to the regime 
where these corrections are small. Therefore the rise of \{tpip) \ for Am >900 MeV should 
not be taken too serious. 

The results of our calculations can be summarized by noting that the structure of the 
pion propagator stays relatively stable for Am ^ A q although sizable changes in m n , f n 
and the quark condensate are found. For larger values of Am instabilities occur, which 
might be related to instabilities of the underlying ground state. 

Unfortunately, for numerical reasons we cannot perform calculations in the exact chiral 
limit, m = 0, where the internal pion propagator has a pole for q 2 = 0. However, 
as mentioned above, one can show analytically for the improved pion propagator that 
£)~ 1 (0) = in this case, independent of Am- Furthermore, it is quite reasonable to expect 
that the slope of D~ x will behave similarly to what we found for mo = 6.1 MeV. This 
implies that above a certain value of Am the residue at q 2 = gets again the "wrong" 
sign, indicating an instability. 

In Fig. |3| we did not change the parameters which were determined in leading order 
by fitting f<?\ and (^) (0) . Of course, if one wants to apply the model to describe 
physical processes a refit of these observables should be performed including the 1/N C - 
corrections. This was partially done in Refs. and || (without taking into account the 
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full momentum dependence of the quark triangles (see Fig. |I] (a)) and the other quark 
loops), where for each choice of Am the new /„. was fitted to the empirical value. Since the 
l/A c -corrections to f n are negative this means that the leading-order term has to be 
larger than the empirical value. As a consequence, the onset of the instabilities is shifted 
to much larger values of Am- In order to estimate the effect, we recall that without refit 
the second pion pole shows up at Am ~ 1250 MeV. At this point we find f n = 35 MeV 
and m n = 216 MeV. A simple way to get the correct pion decay constant f^ xp = 93 MeV 
is to rescale the parameters such that all quantities of dimension energy are multiplied by 
a factor f% xp /fn, which is about 2.7 for the above example. Hence, the point at which the 
second pion pole emerges corresponds to a rescaled cutoff Am — 3300 MeV. Of course, the 
pion mass is rescaled by the same factor and is now much too large (~ 570 MeV). This can 
be easily cured by choosing a smaller current quark mass, which does not affect f n very 
much. Making use of the (almost) linear dependence of on m we obtain m — 1 MeV. 
However, with these parameters we would strongly overestimate the quark condensate. 
Performing a more careful parameter study, it turns out that a simultaneous fit of f n and 
{ipip) is only possible for Am ^ 700 MeV. Therefore the region where the instabilities show 
up in the pion propagator seems to be far away from a realistic parameter set. Ultimately, 
one should of course try to determine also Am itself, e.g. by fitting the decay width of 



vector mesons [11 



4 Conclusions 

We have calculated the pion propagator within an NJL model which was extended to ex- 
plicitly include meson loops in a 1 /Ac-expansion. As already pointed out in Refs. |7J 
and || this requires the introduction of a new cutoff parameter, reflecting the non- 
renormalizability of the NJL model. Similar to these authors an independent meson loop 
cutoff Am has been introduced in addition to the cutoff A q which was used to regularize 
the quark loops. Since the importance of the l/A c -correction terms depends strongly 
on Am, the question of how much the results are altered and in particular whether the 
spontaneous breakdown of chiral symmetry is spoiled by mesonic fluctuations , cannot 
be uniquely answered. Whereas for large values of Am instabilities show up in the pion 
propagator, this is not the case at lower values of A M - First estimates seem to indicate 
that the region of instabilities is quite far away from a realistic parameter set, leaving 
enough room for further applications of the model to physical processes. 

The instabilities in the pion propagator might be a hint for instabilities of the under- 
lying ground state against mesonic fluctuations. However, we should clearly state, that 
this question cannot be assessed within the formalism we have presented here. Similarly 
the stable pion propagator we find at lower values of Am does not prove that the ground 
state is stable. (Of course it should be stable if Am is sufficiently small.) Here further 
careful studies are necessary including a comparison of different regularization schemes. 
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